The Practality of the Hydrodynamic Modes of Diffusion 
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Transport by normal diffusion can be decomposed into the so-called hydrodynamic modes which 
relax exponentially toward the equilibrium state. In chaotic systems with two degrees of freedom, 
the fine scale structure of these hydrodynamic modes is singular and fractal. We characterize them 
by their Hausdorff dimension which is given in terms of Ruelle's topological pressure. For long- 
wavelength modes, we derive a striking relation between the Hausdorff dimension, the diffusion 
coefficient, and the positive Lyapunov exponent of the system. This relation is tested numerically 
on two chaotic systems exhibiting diffusion, both periodic Lorentz gases, one with hard repulsive 
forces, the other with attractive, Yukawa forces. The agreement of the data with the theory is 
excellent. 
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Recent progress has revealed very interesting relation- 
ships between non-equilibrium statistical mechanics and 
and fractal structures that appear in the phase space of 
non-equilibrium, chaotic systems. Theoretical and nu- 
merical work has shown that systems of interacting par- 
ticles have a high sensitivity to initial conditions as char- 
acterized by positive Lyapunov exponents |]l|,|U|^,^,^| . 
As a consequence, the motion of the particles compos- 
ing the system is typically chaotic. Under stationary 
or time-dependent non-equilibrium conditions, this dy- 
namical chaos generates fractal structures in phase space 
^, ppp| P,p^,p^pIp|,p^|,p7t . In this context, Gilbert et 
al. p!8| have very recently shown for a multibaker model 
of diffusion that the entropy production expected from 
non-equilibrium thermodynamics is a consequence of the 
fractal character of the hydrodynamic modes. Moreover, 
the same authors have shown that the diffusion coeffi- 
cient is related to the fractal dimension of these modes 
in a family of multibaker models Jl9| . The important 
question remains whether such results can be extended 
to more realistic systems. 

The purpose of this Letter is to describe the fractal- 
ity of the hydrodynamic modes of diffusion for general 
chaotic systems with two degrees of freedom, including 
Hamiltonian ones. The present approach requires neither 
an escape of particles out of the system as in the escape- 
rate formalism, nor a non-Hamiltonian thermostat as in 
the thermostatted-system approach. Both formalisms 
lead to connections between transport coefficients and 
chaotic properties, but in a restrictive framework, which 
is not the case here. We obtain a formula which gives 
the Hausdorff dimension of the hydrodynamic modes in 
terms of: (a) the dispersion relation of diffusion, and (b) 



the Ruelle topological pressure which characterizes the 
chaotic properties of the dynamics. This new formula 
allows us to prove that, for general chaotic systems with 
two degrees of freedom, the diffusion coefficient is given 
in terms of the Hausdorff dimension of the modes and 
the positive Lyapunov exponent, as conjectured in Ref. 
Jl9| . The results are numerically tested for two examples 
of Hamiltonian Lorentz gases. 

The relaxation rate of a hydrodynamic mode of 
wavenumber k is given by the decay rate of the Van Hove 
intermediate incoherent scattering function pC| as 



s k = lira 7 In (exp [ik ■ (r t - r )]) , 



(1) 



where (•) denotes an average over an ensemble of initial 
conditions and r denotes the position of the tracer par- 
ticle. We assume that the system is mixing so that this 
average weakly converges toward its equilibrium value. 
The diffusion coefficient T> can be obtained by expand- 
ing the relaxation rate (|l]) in powers of the wavenumber 
as 



s k = -Vk 2 + C(k 4 ) 



(2) 



which is known as the dispersion relation of diffusion 
pl[ . In deterministic chaotic systems, the hydrodynamic 
modes of diffusion turn out to be given by a singular 
distribution without density function p^ |. For singular 
distributions, a cumulative function can nevertheless be 
defined by 



F k (0) = lim 



f°dO'exp{ik-{r t (0i)-r Q (6')}} 



J^d0'exp{tk-[r t (O')-r o (9')}} 



(3) 
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where the initial conditions form a one-dimensional line 
parameterized by the angle 8 £ [0, 2tt[. We suppose that 
the system has two degrees of freedom so that each tra- 
jectory has one unstable, one stable, and two neutral 
directions in the four-dimensional phase space. The one- 
dimensional line of initial conditions is assumed to be 
transverse to the unstable directions. The denominator 
in Eq. (||), allows us to define a normalized cumulative 
function such that -Fk(O) = and Fk(27r) = 1. The cumu- 
lative function forms a curve (Re i 7 ^, Im Pk) in the com- 
plex plane C. This curve is invariant under the dynam- 
ics and it gives the fine structure of the hydrodynamic 
mode of relaxation along the given one-dimensional line 
in phase space. 

If we assume that the system is chaotic with the 
Axiom-A properties, we can use a result by Sinai, Bowen 
and Ruelle that averages can be performed in terms of a 
subset of unstable trajectories covering the phase space 
within a certain resolution [^3|,^4) . As the time t and the 
resolution increase the number of these trajectories also 
increases so as to fill densely the phase space in the limit 
t — > oo. The probability weight given to each of these 
trajectories is inversely proportional to their stretching 
factor Aj which characterizes its dynamical instability. 
We notice that the Lyapunov exponent of the trajectory 
is given by A^ = hm{_ >00 (l/t) In |Aj |. Under these cir- 
cumstances, the average of a quantity A can be expressed 
in the limit t -> oo as (A) ~ £\ [A^l^AW, where A^ 
is the quantity evaluated for the j th trajectory. Accord- 
ingly, Eq. (pi) can be transformed into the condition 
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(i)i 



exp(— Ski) exp |ik • 



M „(i) 
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(4) 



Indeed, Eq. (Q) is obtained after taking the logarithm of 
Eq. (^), dividing by t, and taking the limit t —* oo. 

Now, let us suppose that the sum in Eq. (||) is re- 
stricted to the trajectories issued from initial conditions 
in the interval [0, 8}. In this case, we obtain at time t a 
polygonal approximation of the cumulative function Eq. 
(|J) because, the integral L is an average over trajecto- 
ries with initial conditions in [0, 8] and the denominator 
is proportional to the factor cxp(ski). In this respect, we 
notice that Eq. (||) with a complete sum is equivalent to 
the condition Pk(2-7r) = 1. Therefore, we can conclude 
that, at time t, the curve (Re Pk,Im Pk) C C is approx- 
imated by a polygon of sides given by the small complex 
vectors 



AF^ = lA^I" 1 exp(-ski) exp {ik ■ 



-t r o 



}■ 



(5) 



as in the construction of fractal curves of von Koch's type 
p5[ . Each side has the length 



\AF&\ 



exp(— Re Sk t) 



(6) 



so that the polygon can be covered by balls of diameter 
Sj. In the limit t — > oo, this polygon converges to a fractal 
curve characterized by a Hausdorff dimension given by 
J2j zf 11 ~ 1 • Accordingly, the Hausdorff dimension of 
the hydrodynamical mode should satisfy the condition 



exp(— Djj Re Sk t) 



(7) 



On the other hand, Ruelle's topological pressure is de- 



fined in dynamical systems theory by |l3,15 26 



P(J3)= lim - lndAtl 1 "^) 



(8) 



Since there is no escape of trajectories in our case, we 
have that P(l) = 0. The mean positive Lyapunov expo- 
nent of the system is given by A = —dP/d/3\p=i. Equa- 
tion (||) can be transformed in the same way as Eq. (|l|) 
is transformed into Eq. (0), so that 



£|A^V exp[-P(/3) t] ~ f _l 



(9) 



Comparing with Eq. (]?]), we infer that j3 = -Dh where- 
upon 



P(D H ) = D K Re s k 



(10) 



which is the new and central formula of this Letter. This 
formula generalizes Bowen's formula P(Z?h) = for the 
Hausdorff dimension of a fractal invariant set of trajec- 
tories ^?J. In the present case, we are instead dealing 
with a complex invariant curve associated with a relax- 
ation process of rate Sk, which explains the presence of 
a new term in the right-hand side of Eq. (]l(]) . The root 
of Eq. ( |lC| ) gives the Hausdorff dimension Dh (k) of the 
hydrodynamic mode of wavenumber k. 

If the wavenumber k vanishes, the relaxation rate (|l|) 
also vanishes, Sk=o = 0, and Eq. (||) shows that the 
cumulative function becomes -Fk=o(#) = 0/(27r), which 
forms a straight line in the complex plane. In this equi- 
librium limit, Eq. (|ic| ) reduces to P(-Dh) = so that we 
get the dimension Dh = 1, as it should be for a straight 
line. 

For a non-vanishing but small wavenumber, the Haus- 
dorff dimension is expected to deviate from unity. Insert- 
ing Dh = 1+S and the dispersion relation (||) in Eq. (|ic|), 
we can expand both sides in powers of the wavenumber 
by using the aforementioned properties of the topological 
pressure. This straightforward calculation shows that the 
Hausdorff dimension of the hydrodynamic mode is given 
by 



D H (k) = 1 



V 



0(k 4 ) 



(11) 
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as previously obtained for the multibaker models [19|. 
Accordingly, we have proved the conjecture of Ref. [19|, 
that for general chaotic systems with two degrees of free- 
dom, the diffusion coefficient can be expressed in terms 
of the positive Lyapunov exponent A and the Hausdorff 
dimension of the hydrodynamic modes as 



V 



A lim 



k 2 



(12) 



We numerically tested the relation ( |l2| ) for periodic 
Lorentz gases, either with repulsive or with attractive 
scatterers. For the case of repulsive scatterers, we con- 
sidered the periodic Lorentz gas where a point particle 
undergoes elastic collisions on hard disks forming a tri- 
angular lattice. This billiard is known to be fully chaotic 
p8[ . If the horizon seen by the moving particle is finite 
the diffusion coefficient is positive and finite [ ^8|j29[| , and 
the higher-order coefficients such as the super-Burnett 
coefficients are finite ||(]] . We take the radius of the disks 
to be unity, and the distance between the centers of the 
disks equal to d — 2.3, so that the finite-horizon condi- 
tion 2 < d < 4/V3 is satisfied. Numerical computation 
shows that the diffusion coefficient is T> = 0.25 ± 0.01 
and the positive Lyapunov exponent A = 1.76 ± 0.01. 
For this system, the curves (Re F^, Im F^) of the hydro- 
dynamic modes of wavenumbers k x — 0.0, 0.5, and 0.9 
and k y = are depicted in Fig. [p. The initial conditions 
of the point particle are taken as the one-dimensional line 
given in phase space by the positions ro (9) = (cos 9, sin 9) 
and the velocities Vq(8) — (cos 8, sin 8). The Hausdorff 
dimension of the fractal curves of Fig. [I] has been com- 
puted by a box-counting algorithm for different values of 
the wavenumber. These values are plotted versus k x in 
Fig. |^ (filled circles) together with a solid line of slope 
T>/\ = 0.14. The very nice agreement shows that the 
dimension behaves as expected from Eq. (|TT|), which 
provides numerical evidence of the validity of Eq. (|l2|). 




Re F. 



FIG. 1. Periodic Lorentz gas where a point particle of unit 
mass and velocity undergoes elastic collisions on hard disks 
of unit radius forming a triangular lattice with interdisk dis- 
tance d = 2.3: Curves of the cumulative functions of the hy- 
drodynamic modes of wavenumber k x = 0.0, 0.5, and 0.9 with 



k u 



■ 0. Note that the fractality increases with the wavenum- 



ber. The curves are constructed by averaging over 10 points. 




FIG. 2. Hausdorff dimension Dh of the hydrodynamic 
modes versus k 2 = fc 2 . (k y — 0) for both periodic Lorentz 
gases with hard-disk scatterers (filled circles) and Coulomb 
scatterers (open circles). Both solid lines have slopes equal to 
T>/\ for the respective diffusion coefficient T> and Lyapunov 
exponent A of the Lorentz gases. 

For the case of attractive scatterers, we have con- 
sidered the periodic Lorentz gas where an electrically 
charged particle of mass m moves in a square lattice of 
screened Coulomb potentials (also known as Yukawa po- 
tentials) . The Hamiltonian of this system is 



H 



p; 



P 2 y 



2m 



E 



exp(- 



III 



HI 



(13) 



where r = (x, y), 1 = (l x , l y ) with l x , l y = 0, ±1, ±2, ±3, ... 
and || • || is the Euclidean distance. This system was 
proved by Knauf to be fully chaotic and diffusive if 
the energy of the moving particle is positive and large 
enough pl| . We have here taken the parameter values 
m = 1, a = 2, and the energy E = 3. We have numer- 
ically computed the diffusion coefficient V = 2.5 ± 0.1 
and the Lyapunov exponent A = 6.2 ± 0.1. The curves 
(Re F^, Im F^) of the hydrodynamic modes of wavenum- 
bers k x — 0.0, 0.2, and 0.4 and k y = are depicted 
in Fig. |3|. The initial conditions of the point particle 
are taken on a circle centered around the scatterer at 
x = y = 0, with positions i"o(#) = 0.25(cos#, sin#) and 
velocities vo(#) = «o(cos#, sin#) such that E = 3. The 
Hausdorff dimension of these fractal curves are also plot- 
ted in Fig. H (open circles), now together with a solid 
line of slope V/\ = 0.40. Here again, we observe a good 
agreement between both the numerical data and the the- 
oretical prediction of Eq. ([TT|). In this case, the data 
points deviate from the solid line at large values of fc 2 .. 
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Such a behavior is due to the terms of order k 4 and higher 
in Eq. (11). These terms turn out to be more important 



in the case of the Coulomb scatterers than in the case of 
the hard-disk scatterers. Nevertheless, the data points 
converge to the unit value tangentially to the solid line, 
as predicted by Eq. ([ll]). Therefore, we have here also a 
good confirmation of Eq. ( |l2| ) for the Hamiltonian sys- 
tem (ph. 



E 




Re F 
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FIG. 3. Periodic Lorentz gas where a charged particle of 
unit mass and energy E = 3 moves in a square lattice of 
screened Coulomb potentials with a — 2: Curves of the cu- 
mulative functions of the hydrodynamic modes of wavenum- 
ber k x = 0.0, 0.2, and 0.4 with k y = 0. Note that here also 
the fractality increases with the wavenumber. The curves are 
constructed by averaging over 10 6 points. 



In conclusion, we have derived the formula (10) for 
the Hausdorff dimension of the hydrodynamic modes of 
diffusion in chaotic systems with two degrees of free- 
dom. This formula generalizes Bowcn's formula and it 
allows us to establish in more generality a relation be- 
tween the diffusion coefficient, the Lyapunov exponent, 
and the Hausdorff dimension of the hydrodynamic modes 
as conjectured in Ref. |h|. Our results exhibit new rela- 
tions between transport coefficients and chaotic proper- 
ties for systems relaxing to equilibrium, extending those 
previously obtained in the escape-rate formalis m |]7 ,[l3 ,[l4| 
and in the Gaussian thermostat method [p|Js|JTo|! l"j| 17| . 
Therefore, they illustrate the great importance of such 
relationships for both Hamiltonian and non-Hamiltonian 
systems. Extensions of these and related results to fluid 
systems are now under study by the present authors. 
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